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ABSTRACT 

We present analysis of bulk and random gas motions in the intracluster medium using high-resolution Eule- 
rian cosmological simulations of sixteen simulated clusters, including both very relaxed and unrelaxed systems 
and spanning a virial mass range of 5 x 10'-'-2 x lO'^/i"' Mq. We investigate effects of the residual subsonic 
gas motions on the hydrostatic estimates of mass profiles and concentrations of galaxy clusters. In agreement 
with previous studies we find that the gas motions contribute up to « 5%- 15% of the total pressure support 
in relaxed clusters with contribution increasing with cluster-centric radius. The fractional pressure support is 
higher in unrelaxed systems. This contribution would not be accounted for in hydrostatic estimates of the total 
mass profile and would lead to systematic underestimate of mass. We demonstrate that total mass can be recov- 
ered accurately if pressure due to gas motions measured in simulations is explicitly taken into account in the 
equation of hydrostatic equilibrium. Given that the underestimate of mass is increasing at larger radii, where 
gas is less relaxed and contribution of gas motions to pressure is larger, the total density profile derived from 
hydrostatic analysis is more concentrated than the true profile. This may at least partially explain some high 
values of concentrations of clusters estimated from hydrostatic analysis of X-ray data. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 



L INTRODUCTION 



Clusters of galaxies are powerful cosmological probes and 
have the potential to constrain properties of dark energy 
and i ts evolution (e. g., Henrv & Arnaud 1991; Hold er et al 
20011: iMajumdar & Mohr .2003^ ,Hu .2003: .Albrecht et al 



20061 IVikhlinin et alj|2009bl) . Most of the cosmological ap- 



plications using clusters rely on the estimates of their total 
virial mass — a quantity which is difficult to measure ac- 
curately in observations. Traditionally, the masses are esti- 
mated using observable X-ray properties of the intracluster 
plasma using spectroscopi c X-ray temperature or luminos- 
ity (e.g., iRosati et"ani2002L for a review) or velocity disper- 
sion of cluster galaxies. More r ecently, mass estimates via 
strong lensing jS mith et aL 2005), weak lensing (e.g., Dahle' 
I2006t lMahd"^et al., .2008: Zhang et alJ 1 20081). combined 
strong and weak lensiri g (Bradac et al. 2006: iLimousin et alJ 
2p07;'Oguri et al."2009'), and Sunyaev-Zel'dovich effect (e.g., 
Carlstrom et al. 2002, for a recent review) have become avail- 
able. 

One of the most widely used methods for measuring clus- 
ter masses utilizes gas density and temperature profiles of the 
intracluster medium (ICM) derived from X-ray observations 
to solve the equation of hydrostatic equilibrium (HSE) for the 
total mass profile under assumptio ns of spherica l symmetry 
and equilibrium (e.g., Sarazin 1986i:lEvrardetan[i996 ). With 
the advent of high-resolution observations from the Chandra 
and XMM-Newton satellites, which allowed measurements 
of the density and temperature profiles to large radii, it be- 
came possible to use this method without simplifying assump- 
tions of isothermality and /3-model for gas distribution (e.g.. 
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IVikhlinin etan i2006^. 

Recent tests using cosmological simulations of clusters 
have shown that with resolved profiles and with reliable sub- 
traction of the X-ray background, the HSE method can de- 
rive masses of relaxed clusters with an accuracy of better 
than w 10% -20% (Ras ia et alj 120061: iNagai et all [2007bh . 
At the same time, simulations uniformly show the presence 
of ubiquito us subsoni c flows of gas even in very relaxed 
clusters (Evrard '1990'; 'Norman & Brvan"1999'; 'Nagai et al ' 
2003; Rasia etal. 2004; Kav et al. 2004; Faltenbacher et 
2005; Dolagetal. 2005; Rasia et al. 2006; Nagai et al 
2007b: Piffaretti & Vald arnini 20 08; Jeltema et al., ,200: 



llapichino & Niemeveiil2d0 8: Ameg Uo et alj|2009l) 

Such gas motions are thought to be driven by continuing 
accretion of gas onto clusters along filaments, mergers and 
supersonic motions of galaxies through the ICM. Shocks can 
leave behind wakes with sizes comparable to the length scale 
of the cluster Energy of the large-scale eddies in the flow 
can cascade down to smaller scales. Gas motions on smaller 
scales ca n also be d riven directly by motions of groups and 
galaxies (Kim 2007) and by jets and bubbles from the active 
galactic nuclei ( AGN, Churazoy e t al. 2002), although the lat- 
ter may not be as important energetically as mergers in con- 
tributing to gas motions and their contribution may be con- 
fined to the inner regions of clusters. 

Given that only thermal pressure is taken into account in the 
HSE analysis, presence of random gas motions (or any other 
non-thermal pressure component) can contribute to the pres- 
sure support in cl usters an d bias HSE measurements of the 
total mass profiles jEvrard 1990). Analyses of simulated clus- 
ters show that up to w 10% -20% of pressur e support comes 
from subsonic ra ndom gas motions of gas (Rasia et al.ll2004|; 
iFaltenbacher et a l. 2005; Rasia et al.ll2006h . A recent compar- 
ison of mass from w eak len sing and X-ray HSE mass mea- 
surements by Mahda vTet al.l ([2008) shows that the HSE mass 
is biased low by 20% compared to lensing measurement sug- 
gesting that a non-hyd rostatic component in the gas pressure 
gradient (see, however. IZhang et"ani2008h . 
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Random gas motions can alter the structure of the ICM 
through the redistribution of energy from the decay of large- 
scale gas flows with a time scale on the order of a few turnover 
time of major mergers. In particular, it can also be responsi- 
ble for dispersing metals from the ICM core, where the abun- 
dance profile is broader tha n the central galaxy brightness 
profile dRebusco et alj l2005h and more generally in mixing 
gas of differen t metallicity and entropy dWadslev et al.ll2008c 
[Mitchell et al.l l2008). Incomplete thermalization of gas mo- 
tions results in lower ICM temperatures for clusters of a given 
mass, which can contribute to the bias and scatter in the clus- 
ter mass-temperature relation. Viscous dissipation of gas mo- 
tions can lead to secular heating of the ICM. Random gas mo- 
tions can also mainta in and amplify cluste r magnetic fields via 
dynamo processes ( Roettiger et al J [19991: ISubramani an et all 
I2OO6) and cont ribute to the acceleratio n of cosmic rays in the 
ICM (e.g., iBrunetti & Lazarianll2007h . It is therefore impor- 
tant to understand gas flows in the ICM in order to understand 
both thermal and non-thermal processes in the ICM. 

Direct measurements of gas velocities via X-ray spec- 
troscopy (Inogamov & Sunvaev 2003) are challenging with 
the current X-ray instruments, but there are indirect indica- 
tions that residual gas motions ar e indeed present in the I CM 
of observed clusters. In particular. lSchuecker et al.l (|2004|) use 
the Fourier analysis of fluctuations in the projected gas pres- 
sure map of the Coma Cluster measured with XMM-Newton 
and find a Kolmogorov-like spectrum. Given that the prop- 
erties of gas motions and the Reynolds number of the ICM 
in observed clusters are poorly known, cosmological simula- 
tions of cluster formation remain the best tool for studying 
properties of gas motions and evaluating their possible effect 
on observable properties of clusters. 

In this paper we focus on quantifying the contribution of 
random gas motions to pressure support in clusters and the 
corresponding bias in the HSE estimates of total mass pro- 
file using a suite of high-resolution cosmological simulations 
of cluster formation. Our simulations properly treat both col- 
lisionless dynamics of dark matter and stars and gasdynam- 
ics in a self-consistent cosmological setting and capture a 
variety of physical phenomena from the nonlinear collapse 
and merging of dark matter to shock-heating and radiative 
cooling of gas, star formation, chemical enrichment of the 
ICM by supernova and energy feedback. These simulations 
should therefore faithfully capture the dynamical processes 
driving ICM motions and can give us useful insights into its 
expected effects. Although a number of recent studies have 
examin ed random gas rnotions and their effect on mass es- 
timate jRasia et al.l l2004t iKav et alj|2004l: iDolag et alj|200l 
iRasia et al.ll2006l) . most of the studies have used simulations 
with SPH gasdynamics. The magnitude and effects of gas mo- 
tions in such simulati ons depe nds on the specific treatment of 
artificial viscosity (Dolag et al. 2005). Our study, which em- 
ploys simulations with Eulerian gasdynamics with very low 
numerical viscosity, is therefore useful in evaluating possible 
differences between numerical techniques and systematic the- 
oretical uncertainties. 

In agreement with previous studies, we find that gas mo- 
tions contribute up to > 5 - 15% of the pressure support and 
leads to bias in the HSE mass measurement of a similar mag- 
nitude. This pressure contribution and mass bias increase with 
cluster-centric radius and for unrelaxed systems. The increas- 
ing underestimate of mass at larger radii results in derived 
total density profile that is more concentrated than the true 
profile in the simulation. This leads to an overestimate of 



concentrations by 20% when the deri ved profiles are fit with 
the NEW profile ( iNavarro et alj|1996l) . This can at least par- 
tially explain high concentrations found in some recent ob- 
servational studies (Maughan et al. 20071 iBuote e t al. 2007). 
We also demonstrate that total mass and density profiles can 
be recovered accurately if random gas pressure is explicitly 
taken into account. 

The paper is organized as follows. In § |2] we present and 
describe our simulation sample of clusters. In Section 3 we 
analyze the velocity structure of gas and dark matter in the 
simulated clusters, and give the results of the relative fractions 
of random gas pressure and its gradient, and the expected bias 
in hydrostatic cluster mass estimation. In Section 4 we sum- 
marize and discuss our findings. 

2. THE SIMULATIONS 

In this study, we analyze high-resolution cosmological sim- 
ulations of 16 cluster-sized systems in the flat ACDM model: 
= 1 - f^A = 0.3, f^b = 0.04286, h = 0.7 and erg = 0.9, where 
the Hubble constant is defined as 100/ikms"' Mpc"', and 
an (Tg is the power spectrum normalization on an 8/i~' Mpc 
scale. The simulations were done with the Adaptive Refine- 
ment Tree (ART) A^-b ody+gasdynamics code ( lKravtsovll999t : 
iKravtsov et al.ll2002l) . a Eulerian code that uses adaptive re- 
finement in space and time, and (non-adaptive) refinement in 
mass (Klvpinetal. 2001) to reach the high dynamic range 
required to resolve cores of halos formed in self-consistent 
cosmological simulatio ns. The simulation s pre sented here 
are disc ussed in detail in lNagai et al.l (l2007bh and lNagai et alj 
(I2007al) and we refer the reader to these papers for more de- 
tails. Here we summarize the main parameters of the simula- 
tions. 

High-resolution simulations were run using a uniform 128-^ 
grid and 8 levels of mesh refinement in the computational 
boxes of 120/1-1 Mpc for CLlOl-107 and 80/!"' Mpc for 
CL3-24 (see Table [l]). These simulations achieve a dynamic 
range of 32768 and a formal peak resolution of w 3 .66 /z"' kpc 
and 2.44/1"' kpc, corresponding to the actual resolution of 
« 7/1"' kpc and 5/i"' kpc for the 120 and 80/!"' Mpc boxes, 
respectively. Only the region of '--^ 3- lO/i"' Mpc around the 
cluster was adaptively refined, the rest of the volume was fol- 
lowed on the uniform 128^ grid. The mass resolution, mpan, 
corresponds to the effective 512^ particles in the entire box, 
or the Nyquist wavelength of Anj = 0.469/!"' and 0.312/!"' 
comoving Mpc for CLlOl-107 and CL3-24, respectively, or 
0.018/!"' and 0.006/!"' Mpc in the physical units at the ini- 
tial redshift of the simulations. The dark matter (DM) particle 
mass in the region around the cluster was 9.1 x 10** /i"' Mq for 
CLl 1-1 07 and 2.7 X 1 0** /!"' Mq for CL3-24, while oflier re- 
gions were simulated with lower mass resolution. 

To gauge the possible effects of resolution on our results, 
we have used re-simulations of one of the clusters (CL6) with 
varying maximum refinement levels from 6 to 9 (correspond- 
ing to difference in peak resolution of a factor of 8). 

The A^-body+gasdynamics cluster simulations used in this 
analysis include coUisionless dynamics of DM and stars, gas- 
dynamics and several physical processes critical to various as- 
pects of galaxy formation: star formation, metal enrichment 
and thermal feedback due to supernovae Type II and Type la, 
self-consistent advection of metals, metallicity dependent ra- 
diative cooling and UV heating due to cosmological ionizing 
background ( IHaardt & Madaulll996 ). The cooling and heat- 
ing rates take into account Compton heating and cooling of 
plasma, UV heating, and atomic and molecular cooling, and 
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TABLE 1 

Properties of the simulated clusters at z = 0. 



Cluster ID ^soor ''500c ''sooi- relaxed (1) 

(10'-^ A-'Mq) (/r' Mpc) (kms-') unrelaxed (0) 

CLIOI 9.02 1.16 1830 

CL102 5.45 0.98 1547 

CL103 5.70 0.99 1571 

CL104 5.40 0.98 1543 1 

CL105 4.86 0.94 1489 

CL106 3.47 0.84 1332 

CL107 2.57 0.76 1205 

CL3 2.09 0.71 1125 1 

CL5 1.31 0.61 964 1 

CL6 1.68 0.66 1046 

CL7 1.42 0.63 989 1 

CL9 0.83 0.52 826 

CLIO 0.67 0.49 770 1 

CLll 0.90 0.54 847 

CL14 0.77 0.51 806 1 

CL24 035 039 619 

are tabulated for the temperature range 10^ < T < 10"^ K and 
a grid of metallicities, an d UV intensi ties using the Cloudy 
code (ver. 96b4; Ferland etal]|1998h . The Cloudy cooling 
and heating rates take into account metallicity of the gas, 
which is calculated self-consistently in the simulation, so that 
the local cooling rates depend on the local metallicity of the 
gas. Star formation in these simulations was done using 
the observationally-motivated recipe (e.g., Kennicutt 199^: 
P* - Pgas/^*' ^i'^h = 4 X 10^ yrs. The code also accounts for 
the stellar feedback on the surrounding gas, including injec- 
tion of energy and heavy elements (metals) via stellar winds, 
supernovae, and secular mass loss. 

These simulations therefore follow the formation of galaxy 
clusters starting from the well-defined cosmological initial 
conditions and capture the dynamics and properties of the 
ICM in a realistic cosmological context. However, some po- 
tentially relevant physical processes, such as AGN bubbles, 
magnetic field, and cosmic rays, are not included. Therefore, 
the simulated ICM properties are probably not fully realis- 
tic in the innermost cluster regions, where these processes 
are likely important. Moreover, the gas in the simulations 
is treated as an ideal inviscid fluid with a small amount of 
numerical Lapidus viscosity, and it remains unclear to what 
extent the level of motions in the ICM found in the simula- 
tions and discussed below applies to real clusters. Despite 
these limitations, the current simulations reproduce the ob- 
served ICM profiles outside cluster cores dNapai et al.ll2007ah 
and are therefore sufficiently realistic for a purpose of the cur- 
rent study. 

Our simulated sample includes 16 clusters at z = and their 
most massive progenitors at z = 0.6. The properties of simu- 
lated clusters at z = are given in Table [T] The total cluster 
masses are reported at the radius rsooc enclosing overdensities 
with respect to the critical density at the redshift of the output 
(below, we also use a higher overdensity level, 2500). The 
corresponding velocity 

lGM«rsooc) 

V5{)0e = \ (1) 

V '"500f 

is the circular velocity at rsooe- 

^ This term is much smaller than the usual artificial viscosity term em- 
ployed in the SPH simulations and is not important dynamically. 
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Fig. 1. — Profiles of the ICM mean radial and tangential velocity, radial 
velocity and tangential velocity dispersion, the anisotropy parameter (3, and 
ratio of the ID velocity dispersion to the sound speed of gas (Mach num- 
ber). The solid lines represents relaxed clusters and dashed lines represents 
unrelaxed clusters. 



The gas motions and mass measurement biases are expected 
to depend on the dynamical state of clusters. We therefore 
consider dynamically relaxed and non-relaxed clusters sepa- 
rately in this work. Following Nagai et al. (2007b), our re- 
laxed subsample is identified based on the overall structural 
morphology of their Chandra images, mimicking the proce- 
dure used by observers. Specifically, we visually examine 
mock 100 ksec images and identify "relaxed" clusters as those 
with regular X-ray morphology and no secondary maximal 
and minimal deviations from elliptical symmetry. By con- 
trast, "unrelaxed" clusters are those with secondary maxima, 
filamentary X-ray stru ctures, or sig nificant isophotal centroid 
shifts (see Figure 1 of Nagai et al. ('2007b) for the typical ex- 
amples of systems classified as relaxed or unrelaxed). 

In order to assess the effects of gas cooling and star for- 
mation, we also repeated each cluster simulation with only 
the standard gasdynamics for the baryonic component with- 
out radiative cooling and star formation. We will use labels 
'coolingH-SF' (CSF) and 'non-radiative' (NR) to refer to these 
two sets of runs, respectively. Our main analyses are based on 
the CSF runs. We present results of the NR runs only in the 
SectionlHand Fig.l$]and[TOl 

3. RESULTS 

3.1. Velocity Profiles 

Following iFaltenbacher et all (1200 5^ we analyze the veloc- 
ity profiles of ICM for the sample of 16 simulated clusters at 
z = 0. Our definition of the cluster center differs from that 
of IFaltenbacher et al.1 (l2005h : we use the location of the par- 
ticle with the largest local matter density. The mean veloci- 
ties in spherical coordinates {vr,ve,v^) on each radial shell are 
computed as mass-averages of gas velocities randomly sam- 
pled on the shell where we take the rest frame of the shell to 
be the peculiar velocity of the cluster (defined as the mass- 
weighted average of DM particle velocities enclosed within 
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'"500e)- Changing the rest frame to that of the gas shell does not 
introduce significant changes to our analyses. We also rotate 
the gas shell such that the z-axis of the shell is aligned with the 
gas angular momentum vector of the shell. We quantify the 
random gas motion using the trace of the velocity dispersion 
tensor dBinney & Tremaine 2008) : 



(2) 



where ij represent the spatial coordinate indices and the over- 
line denotes mass-average over the gas shell. Throughout this 
paper, we take random gas motions to be synonymous with 
random gas motions. 

While we find no significant variation in the velocity pro- 
files with cluster mass, the gas velocity structure between the 
relaxed sample and the unrelaxed sample are different. In the 
top panels of Figure [T] we show the profiles of mean gas ra- 
dial velocity (vv) and the mean gas rotational (tangential) ve- 
locity V,. For both relaxed and unrelaxed samples, the mean 
radial velocity is negative and small (|v,.| < O.OSvsoot) up to 
rsooc- Beyond this radius, the radial velocities become increas- 
ingly more negative with radius, indicating the radial infall of 
gas into the cluster potential. The mean rotational velocity 
is small 0.2v5ooe at r = O.Srsooe) but a bit larger than the 
mean radial velocity, and decreases with radius for both re- 
laxed and unrelaxed samples. The middle panels of Figure [T] 
show the random motions of the radial and rotational com- 
ponents. Random motions in the unrelaxed sample is in gen- 
eral higher overall than the relaxed sample in both velocity 
components. Radial random motion increases with radius for 
both samples, and the relaxed sample rises more steeply than 
the unrelaxed one. Rotational random motions also increases 
with radius but less rapidly that they are approximately con- 
stant over the radial range shown. Also notice that the random 
motions in both radial and tangential components are about 10 
times higher than their mean counterparts. 

For both samples the internal gas flow remains subsonic 
throughout the clusters, but the Mach number ^ of the 
flow increases with radius. Shown in the bottom-right panel 
of Fig. [T] for unrelaxed clusters the Mach number increases 
from ^ 0.3 at r = 0.4rgooc to 0.5 at rsooc and ~ 0.65 at 
2r5ooc- The Mach number for relaxed clusters is in general 
low compared to the unrelaxed clusters, and it increases from 
X ~ 0.15 at r = 0.3rgooc to 0.4 and 0.55 at r = rsooc and 
2?-5ooc, respectively. 

Some level of anisotropy is expected in reality. It is 
therefore important to examine the radial component and its 
tangential component o f the 3D velocity dispers ion. The 
anisotropy parameter (} (iBinnev & Tremain3l2008h 



/3(r)=l- 



2a2(r)- 



(3) 



provides a useful measure of the relative importance of the 
radial and tangential velocity components. There is a signifi- 
cant difference in gas velocity anisotropy between the relaxed 
and unrelaxed sample. In relaxed clusters, near the core of 
the cluster the gas velocity is slightly more tangential, with 
(3 « -0.8 at r « 0.2r500f . The gas velocity becomes increas- 
ingly radial as one goes further away from the cluster center, 
with (3 rising to nearly zero at rgooc. Unrelaxed clusters, on 

t The Mach number M of the flow is given by the ratio of the one- 
dimensional gas velocity dispersion (ctid = ^ {a^ + a})/'i) to the sound speed 

of gas (Csoiind = V^-Pthcr/Pgas)- 
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Fig. 2. — Top: Ratio of pressure from random gas motions to total pres- 
sure as a function of radius. Relaxed clusters are represented by solid lines 
while unrelaxed clusters are represented by dashed lines. Bottom: Ratio of 
pressure from random gas motions gradient to total pressure gradient for the 
16 clusters. Relaxed clusters are represented by solid lines while unrelaxed 
clusters are represented by dashed lines. 



the other hand, exhibit isotropic gas velocity at small radii 
and slowly becoming more mildly radial as radius increases. 

3.2. Pressure From Random Gas Motions in Clusters 

To gauge the relative importance of pressure support from 
random gas motions in clusters, we measure the ratio of the 
isotropic pressure from random gas motions Prand, given by 



_ 1 



(4) 



to the total pressure, Ptot = Pth + Pmnd, as a function of the 
cluster-centric radius for all 16 clusters, where fth is the ther- 
mal pressure of gas. Note that random gas pressure and its 
gradient is sensitive to small-scale clumps and any pressure 
inhomogeneity, and these sources could potentially bias the 
measurements of the pressure gradient and hence the hydro- 
static mass estimate. To minimize such bias, we remove 
subhaloes with mass greater than 10'^ /i"'M0 and the mass 
within their tidal radius from the calculation. Lowering this 
mass threshold further to 10^ /;~'M0 (which excludes all 
subhaloes) or increasing the radius around subhaloes within 
which material is removed produces no differences in our re- 
sults. In addition, we smooth the pressure and pressure gradi- 
ent profiles using a Savitzky-Golay (SG) filter with a smooth- 
ing scale / = 0.95r where r is the distance from the cluster 
center and second order polynomial interpolation dPress et alj 
|1992|) . The choice of / and the polynomial order was a re- 
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suit of experimentation and allows to reduce the noise in the 
profiles without oversmoothing them. 

In the upper panel of Figure |2] we plot the fraction of the 
pressure from random gas motions to the total pressure as a 
function of radius for each individual cluster For both re- 
laxed and unrelaxed samples, pressure from random gas mo- 
tions fraction increases with radius for r > O.Srgooe- For re- 
laxed clusters, the pressure from random gas motions fraction 
spans a range from 6% to 15% at rgooc- Unrelaxed systems, 
on the other hand, exhibit a greater pressure from random gas 
motions fraction on average, and with larger scatter ranging 
from 9% to 24% at rgooc- (See also FigureE]) 

The stability of the ICM in the cluster potential is provided 
by the gradient of the gas pressure, and it is this gradient that 
is directly relevant to the hydrostatic mass estimates. Assum- 
ing spherical symmetry, in the bottom panel of Figure |2] we 
plot the ratio of the pressure from random gas motions gra- 
dient {dPrind/dr) to that of the total pressure (dPtot/dr) for 
all clusters. For relaxed clusters, the pressure gradient frac- 
tion is relatively constant within O.Srsoor and increases slowly 
with radius. At rsocv, random gas motions contributes to about 
4%- 10% of the total pressure gradient. For unrelaxed clus- 
ters, the pressure gradient is larger, with a range of 3% to 
18%. Differences among relaxed and unrelaxed systems are 
more pronounced in the gradient than the pressure itself. 

As we showed in § 13.11 the random gas motions are gen- 
erally anisotropic. The isotropic pressure from random gas 
motions expressed in terms of radial velocity dispersion and 
the anisotropy parameter is 

(5, 



This means a purely radial pressure from random gas motions 
is underestimated by 2/3/3 if we assume isotropy. At the scale 
of our interest r ~ rsooc where the anisotropy parameter is 
nearly zero for both relaxed and unrelaxed samples, the pres- 
sure fraction profiles remain very similar to Figure|2]when we 
use radial random motions to calculate pressure from random 
gas motions. 

The relative contributions of the pressure from random gas 
motions and its gradient at r^^^ as a function of cluster mass 
Mgoof are shown in Figure |3] which shows that there is no 
statistically significant trend with mass for both. The figure 
illustrates again the fact that in relaxed clusters random gas 
motions has smaller contribution to the total pressure gradient 
compared to the unrelaxed systems. If gas motions is indeed 
stirred up by dynamical interactions between haloes, then ran- 
dom gas motions should have subsided for relaxed systems 
via dissipative relaxation processes. 

3.3. Effect of random gas motions on the hydrostatic mass 
estimates 

In this section we assess the effects of random gas mo- 
tions on the hydrostatic mass estimate of galaxy clusters. We 
treat the random gas motions using approach similar to treat- 
ment of d ynamical coUisionless sys tem using the the Jeans 
Equation ( Binney & Tremain3l2008 !). The mass of the cluster 
within a radius r can be split into separate components corre- 
sponding to different terms of the Jeans equation: 

Mtot(< r) = Mth+Mrand+M„t + Mst,eam+Mc,oss- (6) 

The support from thermal pressure of the ICM (the HSE mass) 
is 



Mth« r) = 
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Fig. 3. — Top: Plot of pressure from random gas motions fraction at 
''sooc versus cluster mass within rsooc- Bottom: Plot effraction of pressure 
from random gas motions gradient at rsooc versus cluster mass within rjooc- 
Relaxed clusters are represented by filled circles while unrelaxed clusters are 
represented by empty circles. 



the support from random gas motions is 



M,-and(< r) = 



Gpgas V 9r 



(2^2 -af); (8) 



the support from non-random gas rotation is 

-2 

rvf 

Mrot« r) = 
the "streaming" term is given by 

r f _ dvr _ dvr Va, dvr 
G \ or oO sm 9 dm 



(9) 



(10) 



and the "cross" term 



1 d(Ps,,<7^g) ^ 1 d{p^,^cr;^) 
r do rsinO d<f> 



(11) 



Gpgas dr 



(7) 



The last two mass terms are much smaller than the other terms 
and can in principle be neglected, but we include them in 
the actual computation for completeness. Note that spherical 
symmetry for the gravitational potential and steady state are 
assumed in deriving the expression for mass from the Jeans 
equation, and all the physical quantities at a given radius are 
averages over a radial shell. 

Figure H] compares the true and estimated averaged mass 
profiles normalized by M500C for relaxed and unrelaxed clus- 
ters. The comparison shows that the HSE mass profile Mther 
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Fig. 4. — Averaged mass profiles M(< r) normaHzed by M500C. for relaxed 
clusters (top) and unrelaxed clusters (bottom). The solid line shows the actual 
mass profile from simulation, the long dashed line shows the mass profile 
from hydrostatic equilibrium including random gas and thermal pressure, and 
the short dashed line shows the mass profile from hydrostatic equilibrium 
including thermal pressure only. Hashed region shows the \-a eiTor of the 
mean. 



(short-dashed line) underestimates the true mass profile mea- 
sured directly in simulations (solid line). The mass profile 
from the full mass expression with both thermal and gas mo- 
tion support (long-dashed line) recovers the true cluster mass 
profiles quite well for relaxed clusters. For unrelaxed clusters, 
the non-thermal terms recover the true mass well within rsooc, 
but overestimates the true mass beyond rsooc, because the as- 
sumptions of spherical symmetry does not hold true at large 
radii in most unrelaxed clusters. 

Figure |5] shows the relative contribution of mass terms 
Mth, M,and and M^oi- For both relaxed and unrelaxed sam- 
ples, the thermal component dominates the pressure support 
(~ 80-90% at rsooc)- The support due to gas motions in the 
unrelaxed sample is higher than that of the relaxed one. For 
both samples, the rotational term is greater than the random 
gas term in the inner region (r < 0.6r50()c for relaxed sample 
and r < 0.2r5ooc for unrelaxed sample), and decreases with 
radius. The random gas term on the other hand increases 
with radius, and accounts for most of the non-thermal gas 
su pport beyon d rsoof . This generally agrees with the findings 
of iFang et all 12009.) which use the same set of clusters, de- 
spite our assumption of spherical symmetry which they do not 
adopt. 

The biases in the hydrostatic mass estimate at rsoor as a 
function of cluster mass are shown in Figure |6] The bias 
(AM/M) is defined as a fractional difference between the es- 
timated mass (Mest) and the true mass (Mtrue) measured di- 
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Fig. 5. — Averaged mass fraction of various pressure contributions for re- 
laxed clusters (top) and unrelaxed clusters (bottom). The solid line shows the 
contribution from thermal pressure, the long dashed line shows the contribu- 
tion from pressure from random gas motions, and the short dashed line shows 
the contribution from rotational support. Hashed region shows the l-cr error 
of the mean. 



rectly in the simulations. Following iNagai et alJ (|2007b'), the 
biases are computed in two ways for each cluster First, both 
Mest and Mtme are measured in the same physical region en- 
closed within "true" radii (rn-ue) measured directly in simula- 
tions. The bias computed this way is indicated with circles. 
In practice, additional biases in the estimated masses (Mest) 
could arise from a bias in the estimation of a cluster virial ra- 
dius. Therefore, we also compute the bias in the estimated 
hydrostatic mass, Mest(< /"est), enclosed within the estimated 
virial radius, rest, indicated with triangles. Both bias values 
are listed in Table |2] at z = and z = 0.6. For relaxed clus- 
ters, we find that the hydrostatic mass at rsooc is biased low 
by 7% and 11% at z = and 0.6, respectively. The bias is 
smaller if the hydrostatic mass is measured in the inner region 
of clusters (e.g., r25()0f). Both bias and scatter are larger for 
unrelaxed systems. We find no apparent trend of the bias with 
cluster mass. 

Note that the biases in hydrostatic mass repo rted in this 
work are smaller than the values reported in Nag ai et aP 
(2007b) based on analyses of mock Chandra observations. 
Although these values are consistent within 1-sigma, we find 
systematic deviations in the mean. For the relaxed clusters 
at z = 0, the bi as at r = rsnor is 7% which is smaller than 
13% reported in lNagai et all (l2007bh . Figure [8] illustrates that 
the distribution of the HSE mass in l^agai et al. (2007b) is 
skewed towards more negative bias relative to the results of 
the present work. The offset is smaller in the inner regions 
(^ 1% at r = r25()0c). The main difference is that biases re- 
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TABLE 2 

Bias in the hydrostatic mass estimate of clusters 







Sample 


Bias ± Error (la) * 


redshift 


radial range 


(number of clusters) 


AM«rtrue)/M 


AM«re,,c)/M 






all (16) 


-0.077±0.022 


-0.113±0.033 


z = 


< '"2500c 


relaxed (6) 


-0.061±0.023 


-0.08()±0.031 






unrelaxed (10) 


-0.086±0.032 


-0.133±0.049 






all (16) 


-0.094±0.040 


-0.131±0.062 


z = 


< ''SOOc 


relaxed (6) 


-0.074±0.022 


-0.112±0.038 






unrelaxed (10) 


-0.106±0.064 


-0.142±0.098 






all (16) 


-0.066±0.023 


-0.074±0.041 


z = 0.6 


< '"2500c 


relaxed (6) 


-0.042±0.027 


-0.044±0.052 






unrelaxed (10) 


-0.081 ±0.034 


-0.093±0.059 






all (16) 


-0.116±0.029 


-0.200±0.045 


z = 0.6 


< ''SOOc 


relaxed (6) 


-0.111±0.027 


-0.154±0.039 






unrelaxed (10) 


-0.119±0.044 


-0.227±0.068 
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^The scatter can be obtained from multiplying the eiTor by ' 
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Fig. 6. — Fractional differences between the true mass and the HSE es- 
timated mass, AM/M = (Mest — Mtiuc) / l^tme, as a function of cluster mass 
Msooc- The circles and triangles show the hydrostatic mass evaluated at the 
true and estimated fsooo respectively. The solid and open symbols indicate 
relaxed and unrelaxed clusters. 



ported in lNagai et alJ (l2007bh were determined in the analyses 
of the mock Chandra data using analysis procedures similar to 
those used in observations, while the present analysis uses the 
3D profiles of gas density and temperature measured directly 
from raw simulation data. We find that the difference in the 
inferred bias is primarily due to differences in the tempera- 
tu re profiles sl ope derived from the mock Chandra analyses 
in iNagai et aTl (|2007b) and the 3D mass-weighted tempera- 
ture profiles directly measured from simulations. 

3.4. Effect on Concentration Parameters of Galaxy Clusters 

As seen from Figure |4] the bias in the HSE mass increases 
with radius, indicating that the contribution of the random gas 
pressure increases toward cluster outskirts. This results in 
the total density profile that is steeper than the true density 
profile and hence an overestimate of concentration parameter 
cgoo = r^QQc/rs if the derived density profiles of clusters are 





AM,ooc/M 



0.2 



500c 



Fig. 7. — Distributions of the HSE mass bias, AM/M = (Mcst - 
Mt,uc)/Mt,uc, for a sample of 6 relaxed clusters at z=0. The solid line shows 
the distribution where the HSE mass is calculated directly from the gas ther- 
mal pressure and density profiles. The shaded area shows the distribution 
where the HSE mass i s derive d from mock Chandra analyses of simulated 
clusters in lNagai et alj <2007bD . 



fitted with the NFW profile (i Navarro et allll996h . Figure [8] 
shows the bias in the concentration parameters for our sim- 
ulated clusters. The concentrations were estimated by fitting 
the NFW profile to the derived total density profile over ra- 
dial range 0.1 < r/rsaoc < 1-0. We find that the estimated 
concentration parameter based on the HSE assumption is bi- 
ased high on average by about (23 ± 10)% for relaxed sys- 
tems and (47 ±21)% for unrelaxed systems. The scatter is 
(23 ± 5)% and (41 ± 12)% for relaxed and unrelaxed systems 
respectively. 

Given that the bias in the total mass profile we find from the 
idealized analysis of the 3D cluster profiles is lower than the 
bias derived from the analyses of mock Chandra data from the 
same simulation (see discussion in the previous section), we 
can expect that the bias in concentrations may also be some- 
what larger for the mock data analysis. Therefore we have 
fitted NFW concentrations to the total density profiles in the 
same radial range as above derived from mock X-ray data for 
relaxed clusters used in Nagai et al. (2007b). The concentra- 
tion parameters in this case tend to be overestimated by about 
a factor of two. One particular problem we encountered with 
fits to some of the systems was that errors in estimates of con- 
centrations can become catastrophic if the NFW scale radius 
of the fit was smaller than the minimum radius of the radial 
range used in the fit (which corresponds to concentration val- 
ues of > 10 for our choice of the radial range). 

The upward bias of the concentrations estimated from the 
HSE-derived density profiles due to gas motions can help ex- 
plain the surprisingly high values of conc entrations me asured 
in some of the re cent X-ray analyses (Maughan et al. 20071; 
iBuote etalj|2007h . In particular, Buote et al-: (2007) found 
that for a given mass concentrations derived for groups and 

* The lower limit on the radius is chosen not to include the region of the 
profile significantly affected by the dark matter contraction in response to 
baryon condensation I Gnedin et al. 2004), while the outer radius corresponds 
to the largest radius to which the density profile is determined reliably from 
current X-ray data. 
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Fig. 8. — Fractional difference of concentration C500 derived from the fit to 
the true total density profiles of clusters and to the HSE-derived density pro- 
files in the simulations with cooling for relaxed (filled circles) and unrelaxed 
(open circles) clusters. 

clusters from HSE X-ray analysis are somewhat higher than 
expected in the standard ACDM model with erg « 0.7-0.8. 
The bias due to gas motions can expl ain part of this differ- 
ence.^ As noted by iFang et aTl (|2009|) . the actual values of 
concentration bias may be affected by the overcooling prob- 
lem of cluster simulations. However, the existence and sign 
of the bias are generic. In particular, as we show in the next 
section, qualitatively similar bias is obtained from the simula- 
tions of our cluster sample in non-radiative regime. 

4. EFFECTS OF COOLING ON THE MASS AND CONCENTRATION 

BIASES 

Given that the cosmological simulations suffer from the 
"overcooling problem" - a too high a fraction of cooled, con- 
densed gas produc ed in simulations com pared to observations 
(e.g., see Fig. 2 in lKravtsov et al.ir2009h - it is reasonable to 
ask whether the re sults we presented in the previous sections 
are affected by it. Fang et al. (2009) have independently an- 
alyzed the same set of simulations we use here and argued 
that much of the non-thermal pressure support in the simu- 
lations with cooling is due to rotation of gas in the inner re- 
gions (r < O.Srvir), which should not be trusted because ro- 
tation is overestimated due to the overestimated amount of 
cooling gas. 

Figure|9]shows the mass profiles for relaxed clusters in non- 
radiative re-simulations of the 16 clusters used in our analy- 
sis and contributions to the total mass from thermal pressure 
and rotation and random gas gas motions. The figure shows 
that while the relative contribution of rotational velocity in 
these simulations is lower in non-radiative simulations, the 
total contribution of gas motions to the HSE reconstruction 
of the mass profile is similar to that of the simulations with 
cooling. Given similar degree of gas motions, albeit with a 
different contributions from rotation and random components 
in the inner regions, the HSE-derived mass profile is conse- 
quently also biased low for non-radiative simulations and the 

' Simulations also predict that condensation of baiyons and their conver- 
sion into stars should also increase concentrations compared to the purely 
dissipationless simulations (iRudd et al.„2008i) . which are used to make pre- 
dictions for the concentration-mass relation in a given cosmology. 
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Fig. 9. — Top panel: the mass profiles (top) for relaxed clusters in the 
non-radiative re-simulations of clusters used in our analysis; bottom panel: 
contribution of thermal pressure, rotational and random gas motions to the 
total massprofile for the clusters in the upper panel. The notation is the same 
as Figure|4|and Figure[5] 
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Fig. 10. — Fractional difference of concentrations C5oor derived from the 
true total density profile and the HSE-derived density profiles for the non- 
radiative re-simulations of our clusters. The filled circles represent relaxed 
clusters and the open circles represent unrelaxed clusters. As in the runs with 
cooling, the concentrations of the relaxed clusters tend to be biased high due 
to the bias in total mass, which increases with cluster-centric radius. 
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mass bias at rsooe is similar to that in the simulations with 
cooling. 

In addition, the figure shows that the bias in the mass es- 
timate also increases with cluster-centric radius in this case, 
which means that we can expect bias in the HSE-derived con- 
centrations for the total density profile. Indeed, Figure [TO] 
shows that the concentration bias is qualitatively similar in 
non-radiative simulations to that in runs with cooling: five 
out of six relaxed clusters have cf^, biased high on average 
by w 10-20% with concentration catastrophically overesti- 
mated in one of the clusters (compare this figure to Fig. [8]l. 
The results presented in this section clearly show that conclu- 
sions of this study with regard to the bias in the HSE-derived 
total mass profile and NFW concentrations are general and 
are not determined by the unrealistic cooling of the gas in the 
dissipative simulations. 

5. DISCUSSION AND CONCLUSIONS 

We presented analysis of gas motions in the ICM using 
high-resolution adaptive mesh refinement cosmological sim- 
ulations of a sample of sixteen galaxy clusters spanning the 
virial mass range of 5 x 10'^-2 x lO'^/z"' Mq. Our study fo- 
cuses on the effects of residual gas motions on the estimates 
of the total mass profiles of clusters from the hydrostatic equi- 
librium analysis. We analyze systems that appear morpholog- 
ically relaxed and unrelaxed in mock Chandra X-ray images 
separately to study the effects of dynamical state of clusters 

on the resulting HSE mass profile. 

In broad agreement with previous studies ([EvrardI 1990 ; 
Frenket al. 1999; Rasia et al. 2004; Faltenbacher et al 2005'; 
Rasia et al. 2006; Nagai et al. 2007b; Piffaretti & Valdarnini 
20081; IJeltema et all 12008; la pichino & Niemever 2008D, we 
find that gas motions contribute up to « 5% - 15% of the pres- 
sure support in relaxed clusters, which leads to underestimate 
of the total virial mass in the HSE analysis accounting only 
for the thermal pressure. We have used re-simulations of one 
of the clusters with different maximum levels of refinement 
and find that the magnitude of the pressure support due to gas 
motions at the considered radii has converged. 

On average, the hydrostatic cluster mass estimate is biased 
low by about 6 ± 2% at r = rasooe and 8 ± 2% at r = rsooc for 
relaxed systems, while the biases in unrelaxed clusters are 
about 9 ±3% and 11 ±6% at these radii, respectively. We 
have tested that our results are not affected by resolution by 
analyzing re-simulations of one of the clusters with different 
spatial resolution. The magnitude of the effect is consistent 
with the recent observational evidence of a similar bias based 
on comparison of the HSE and weak lensing derived masses 
dMahdavi et al.i r2008). 

We observe that the average Mach number of the ICM gas 
motions is rather small A4 ^ 0.4 for unrelaxed clusters, which 
implies that gas motions are generally subsonic and thus do 
not dissipate via shocks. This means that it may take sub- 
stantial amount of time to thermalize these motions fully, and 
that the time will depend largely on the physical viscosity 
of the gas. Incomplete thermalization in unrelaxed clusters 
leads to a lower measured X-ray temperature Tx that might 
bias t he M-Tx scaling relation with a higher normalization 
(e.g., iMathiesen & Evrardll2001h . On the other hand, under- 
estimating the cluster mass from the hydrostatic equilibrium 
analysis would bias the M-Tx relation to lower masses at a 
given temperature. This issue can potentially be resolved by 
comparing the HSE-derived masses to cluster masses derived 
from weak lensing surveys (iMahdavi et al.ll2008l;IZhang et all 



l2008H\^khiinin et al.ll2009ah . 

The bias in the total mass profile increases with increasing 
cluster-centric radius radius which results in the HSE-derived 
total density profile that is more concentrated that the true pro- 
file. The best fit NFW concentrations fit to the HSE-derived 
profile therefore results in concentration values generally bi- 
ased high compared to the corresponding fits to the true den- 
sity profiles. In particular, we find that the concentration pa- 
rameters, cgoo = rsoQc/rs, based on the HSE mass profile, is 
biased high on average by about ^ 24%. This bias would 
have to be taken into account for observational studies of 
the conce ntration-mass (c-M) relation based on th e HSE as- 
sumption (Mau ghan et al.ll2007t fBuote et al.ll2007h . We show 
that a similar bias exist in re-simulations of the same clusters 
in the non-radiative regime. Our conclusions therefore are 
general and are not due to the overcooling in th e dissipative 
simulations, as was claimed bv lFang et al.l (l2009l) . 

We demonstrate that the mass profile within rsooc can be 
recovered well if the pressure support due to gas motions is 
explicitly taken into account. On the other hand, the recovery 
is not as good beyond rsoof, with the deviation from the true 
mass profiles depending on the cluster dynamical state and 
configuration of mergers. There are some uncertainties in our 
analysis due to limitations of our simulations. Although we 
show that adding gas motion support recovers the mass of a 
relaxed cluster accurately, we have not considered other po- 
tential sources of non-thermal pressure support, such as cos- 
mic rays and magnetic fields, that can also bias the hydrostatic 
mass measurement. These effects have to be considered sepa- 
rately. Note that the effect of gas motions generally increases 
with increasing radius simply because clusters are less relaxed 
at larger radii. The effects of other sources of non-thermal 
pressure support can have different radial dependence and fu- 
ture observations may potentially exploit such differences to 
differentiate among them. 

In addition, real gas motions in the ICM plasma may also 
be different from the motions in our simulated clusters due 
to a different physical viscosity of the intracluster gas com- 
pared to simulations. Given that higher viscosity would tend 
to damp the gas motions faster, the estimates of their ef- 
fects presented here provide an upper limit to the random gas 
pressure of the ICM. Simulations exploring reasonable sce- 
narios for magneto-hydrodynamics effects, effects of AGN 
feedback, and cosmic rays injection and evolution should pro- 
vide some guidance as to the radial dependence of contribu- 
tions from magnetic fields and co smic rays, as non-thermal 
sources of pres sure s upport (e.g., ISubramanian et al.l 120061 ; 
iPfrommer et"ani2007t Isijacki et al.ll2008h . 

Interestingly, we do not detect a strong dependence of the 
mass bias with cluster mass, which means that the bias would 
affect the normalization of the mass observable relations but 
not its slope. However, the sample of simulated clusters is 
rather small. iJeltema et al.l (l2008i) use a considerably larger 
sample of simulated clusters and do detect a we ak trend of the 
HSE mass bias with cluster mass. In addition. IJeltema et al.l 
(2008) show that it is important to use quantitative measures 
of the dynamical state of the clusters rather than the visual 
classification of the mock X-ray images, as was done in our 
analysis (note, however, that we follow the standard practice 
of such classifi cation in obs ervational X-ray analyses). In 
agreement with IJeltema et al.l (12008). our own analysis shows 
that often clusters that appear morphologically relaxed along 
one projection, would not be classified as relaxed along other 
projections (the reverse is also true, as clusters that are in- 
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trinsically relaxed may seem unrelaxed in projection due to 
projecting groups and clusters along the line of sight). 

Observationally, little is known about gas motions in clus- 
ters except that in many systems there are indirect indications 
of bulk gas motions associated with merg ers and motions 
of the central cluster galaxies (see Markev tch & VikhlininI 
l2007l for a review). ISchuecker et al.l (|2004|) presented evi- 
dence for possible random gas motions in the Coma cluster 
via fluctuations in the pressure map. It is not yet clear, how- 
ever, whether this interpretation is unique and whether the re- 
sults are generic for other clusters. The most direct way of 
measuring gas motions would be to measure broadening of 
line p rofiles of heav y ions in radio (ISvunyaev & Churazovl 
119841) or in X-rays (llnogamov & SunvaevI l2003b" In radio 
wavelengths this requires very deep observations to detect the 
lines in emission or sensitive observations against bright ra- 
dio sources in clusters to detect the lines in absorption (D. 
Marrone, priv. communication). Observations of the lines in 
X-rays require sensitive high resolution X-ray spectrometer. 



which have not yet been available, but may become possi- 
ble with the launch of the ASTRO-H satellited A combina- 
tion of constraints from mass measurements via gravitational 
lensing, possible future direct measurements of gas motions, 
and improved modeling of the ICM physics is thus the way to 
make progress in our understanding of the mass measurement 
biases discussed in this paper. 
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